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Abstract: In the first part of this work, we consider second order supersymmetric 
differential operators in the semiclassical limit, including the Kramers-Fokker-Planck 
operator, such that the exponent of the associated Maxwellian is a Morse function 
with two local minima and one saddle point. Under suitable additional assumptions 
of dynamical nature, we establish the long time convergence to the equilibrium for 
the associated heat semigroup, with the rate given by the first non-vanishing, expo- 
nentially small, eigenvalue. In the second part of the paper, we consider the case 
when the function has precisely one local minimum and one saddle point. We also 
discuss further examples of supersymmetric operators, including the Witten Lapla- 
cian and the infinitesimal generator for the time evolution of a chain of classical 
anharmonic oscillators. 
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1 Introduction and statement of the main result 

The principal purpose of the present paper is to apply the spectral results of [12] to 
obtain a precise information concerning the large time behavior of the heat semigroup 
generated by the semiclassical Kramers-Fokker-Planck operator 

P = y.hd x -V>(x).hd y + l(-hd y + y).(hd y + y), x,y £R n , 7 > 0. (1.1) 

In fact, as in [12], our main result will be valid for a large class of supersymmetric 
second order differential operators, including (11.11) . Physically, the semiclassical limit 
h — > in (11.11) corresponds to the regime of low temperatures. Recall that by 
supersymmetry, we mean the fact that the Kramers-Fokker-Planck operator P can 
be viewed as a Witten Laplacian in degree associated to a certain non-semidefinite 
scalar product in the spaces of differential forms. These relations with the Witten 
complex [22] were exhibited in the works of J. Tailleur, S. Tanase-Nicola, J. Kurchan 
[21] and J. M. Bismut p], using respectively the languages of supersymmetry and 
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differential forms (See also [8] for a quick introduction to the differential form version 
of Bismut and [2]). 

The paper [12], which is a natural continuation of [TTJ, analyzed resolvent esti- 
mates and the low lying eigenvalues of P, assuming that the potential V in (II. ip is 
a smooth real valued Morse function on R n , such that 

d a V = 0(l), \a\>2, (1.2) 

and with 

|W| > 1/C, for \x\ > C > 0. (1.3) 

Assuming furthermore that V has precisely three critical points: two local minima, 
x±i, and one critical point Xq of index 1, it was established in [12] that for C > large 
enough and h > sufficiently small, the operator P has precisely two eigenvalues in 
the disc D(0, h/C) = {z 6 C; \z\ < ^}, /i and /ii, such that /x = and /xi is real 
and of the form 

/it = h ( ai {h)e~ 2Sl/h + a^{h)e- 2S - l/h ) , Sj = V(x ) - V{x 5 ). (1.4) 
Here cij are real with 

(ij(h) ~ a^o + hcij i + . . . cijfl > 0. (1.5) 

Notice that the eigenfunction corresponding to the eigenvalue /i = is the Maxwel- 
lian 

exp (-<t>/h) E L 2 (Rll), <Kx,y) = £ + V{x). (1.6) 

In the case when V — > +oo as x — > oo and V has precisely one local minimum, 
it follows from the results of [11] that in a disc D(0,Ch), C > 1, apart from the 
eigenvalue /xq = 0, the real part of the other eigenvalues is > h/C. In this case, 
precise results describing the behavior of the semigroup exp (—tP/h) for large t, were 
obtained in [llj — the rate of the return to equilibrium, given by the projection onto 
the ground state (11 .6p . is then of the order of magnitude 1. In this work, we shall 
complement this study by analyzing the question of a return to equilibrium in the 
presence of exponentially small eigenvalues, due to the tunneling between the local 
minima. 

The following is the main result of this work, specialized to the case of the 
Kramers- Fokker-Planck operator (11.11) . Here we shall also write P for the m-accretive 
realization of the operator (II. ip — see also section 2 and section 3 in [T2] . 
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Theorem 1.1 Assume that V in (11.11) is a C°° real valued Morse function satisfying 
(11.21) and (11.31) and having precisely 3 critical points: 2 local minima, x±\, and a 
critical point of index 1, so that the disc D(0, h/C) for C > large enough, contains 
precisely 2 eigenvalues of P, fio = and \i\ given in (11.41) . Let IL, be the spectral 
projection associated with the eigenvalue fij, j = 0, 1. Then we have 

Tlj = 0(1), h 0. (1.7) 

We have furthermore, uniformly as t > an<i h — > 0, 

e -^ = n o + e-^ l/h n 1 + 0(l)e- t/c , C>0, zn C{L\L 2 ). (1.8) 

The structure of the article is as follows: In section 2, relying upon the results 
of [12] , [11] , we establish the basic a priori coercivity estimate for the operator P in 
a suitable exponentially weighted space, introduced in [12]. In this sense it can be 
interpreted as an hypocoercive estimate (see e.g. [13], [23]). In section 3 it is then 
quite straightforward to prove Theorem 1.1 in its general form, by combining the 
results of section 2 together with the analysis of [12J. 

The second part of the paper, consisting of sections 4-6 is of a somewhat different 
nature, complementing and extending the previous analysis. In section 4, we study 
the case left out in [12], when the potential V in (II. ip has precisely two critical 
points: one local minimum and a critical point of index one. In this case, is not 
an eigenvalue of P, and the large time behavior of the heat semigroup is governed 
entirely by the first, exponentially small, eigenvalue. In section 5, we give some 
examples and describe a probabilistic framework in which the Witten Laplacian and 
the Kramers- Fokker-Planck operator both arise naturally. Finally, in section 6, we 
describe another example of a supersymmetric operator, for which the question of 
a convergence to equilibrium is of interest, namely a chain of classical interacting 
anharmonic oscillators, coupled to a heat bath. We show how to adapt the analysis 
of [12] to cover also this case. 

Acknowledgments: A part of this work was carried out in May of 2007, when the 
second author was visiting Universite Paris 13 being on leave from UCLA. It is a 
great pleasure for him to thank its Departement de Mathematiques, and in particular 
Alain Grigis, for the extraordinary hospitality and excellent working conditions. His 
best thanks are also due to the Centre de Mathematiques of Ecole Polytechnique for 
the generous hospitality in June of 2007. The partial support of his research by the 
National Science Foundation under grant DMS-0653275 and by an Alfred P. Sloan 
Research Fellowship is also gratefully acknowledged. 
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2 An (hypo-) coercive estimate 



The purpose of this section is to establish an a priori estimate for P, instrumental 
in proving Theorem 1.1. This estimate will imply the exponential decay for the heat 
semigroup exp (—tP/h) in L 2 , when restricted to the kernel of the spectral projection 
corresponding to the eigenvalues /Iq — and /ii in (11.4ft . When doing so, as in [12J, 
rather than working directly with ( 11. lft . we shall consider a broader class of scalar 
real second order non-elliptic non-self adjoint operators on R™. For completeness, we 
shall now recall, following [12], the definition and the main assumptions concerning 
this class. 

Let us consider 

n 1 n 

P = ^ hD xj ° b jA x ) ° hD *k + 2 { c j( x ) hd x 3 + hd xj ° c j( x )) + Po( x ) (2- 1 ) 

j,k=l j=l 

=: P 2 + iP 1 +P , D x . = 

Here the coefficients bj t k, Cj, po all belong to C°°(R n ; R), with 6^ = bkj- Associated 
to P in (12.11) is the symbol in the semiclassical sense, 

P(x,0 =P2(x,£) +ipi(x,£) +Po(x), (2.2) 

n n 

P2(x,0 = ^2 hk( x )€jtk, Pi(x,0 = ^2cj(x)£j, (2.3) 

j,k=l 3=1 

so that pj(x, £) is a real- valued polynomial in £, positively homogeneous of degree j, 
< j < 2. We may notice that p(x, £) coincides with the Weyl symbol of P modulo 
0(h 2 ), locally uniformly. 

As in [12], we shall assume that 



p 2 (x,O>0, Po(x)>0. (2.4) 

Furthermore, we shall impose the following growth conditions, 

8 a x b hk {x) = 0(1), \a\ > 0, (2.5) 

d%Cj{x) = Oil), \a\ > 1, (2.6) 

d: Po (x) = 0(1), \a\ > 2. (2.7) 

From section 3 in [12], we recall that under these assumptions, the graph closure of 

P : 5(R") -> S(R n ), still denoted by P and such that ReP > 0, coincides with the 
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maximal closed realization of P, with the domain given by T>(P) = {u G L 2 ; Pu G 
L 2 }. In particular, this shows that the operator P is m-accretive, and hence, the 
contraction semigroup 

e- tp/h : L 2 -> L 2 , t>0 (2.8) 

is well-defined. 

We proceed next to recall the additional assumptions of a dynamical nature, 
introduced in section 4 of [T2]. Let 

n 

z/(x, 9a,) = ^ Cj(x)d xp (2.9) 

and recall the Hypothesis 4.1 of [T2] : 

The set {s G M n ; p (a?) = 0, i/(a;, d x ) = 0} is finite = {x u . . . x N }. (2.10) 
With pj = (xj, 0), 1 < j < N, we define the critical set 

c = {ft,..Mc» 2 " (2.11) 

The coefficients p , pi, p 2 in (12. 2p all vanish to the second order at each pj G C. As 
in [12], we define 

p(x,0=Po(x) + (0~ 2 P2(x,0, (2-12) 
and consider the time average 

I fTo/2 

— / poexp(tH Pl )dt, T >0. (2.13) 

J J-T /2 

We shall assume that the Hypothesis 4.3 of [12] holds true: 

For T > Ofixed, we have near eachp^, (p) To (p) ~ \p — pj\ 2 , (2-14) 
and in any set of the form \x\ < C, dist(p, C) > 1/C, we have 

(p) To {p)>J— C(C)>0. (2.15) 
C{C) 

We also need an additional dynamical hypothesis near oo in M n , 

V neighborhood U of ir x C, and V x G R n \ U, 3C > 0, 



/To 
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Under the assumptions above, the paper [12] defines an auxiliary real valued weight 
function i/j £ (x,£) = 0(e) on T*R n , e > 0, such that 

together with an associated canonical transformation 

k(6) :R 2n -*A s :={(x, £) +i6H^(x,£)) (x,£) el 2 "}, < 6 < 1, (2.17) 
for which 

k(6)(x,0 = (x,^+idH^(x,0+O(e 1 / 2 5 2 ). (2.18) 

Here we let Hf denote the Hamilton field fL(x,£) • d x — f' x {x, £) • <9g of a C^-function 
f(x,£). We refer to section 4 of [12] for the details of the construction of ipe an d Kg. 
Here we shall merely recall that 5 > fixed in ( 12.181) should be small enough, and 
e = Ah, with A arbitrarily large but fixed. 

Associated to k(8) in (I2.17P there is an elliptic Fourier integral operator with a 
complex phase As >e , constructed in Section 5 of [12], such that 

A s , £ :S^S (2.19) 

continuously, and 

As, £ = O a {1) : L 2 -> L 2 . (2.20) 

Moreover, it is proved in [12] that Ag j£ is invertible, when e/h^> 1, with the inverse 
Ajl also enjoying the mapping properties (I2.19p . (12.201) . 

When B > A fixed is to be chosen, and A ^> B is large enough, as in sections 
6,7 of [12], we shall consider the conjugated operator 

= e = Ah, (2.21) 

acting on L 2 (IR n ). It was then proved in [12] that the real part of the symbol of P s,£ 
is > § - Ch outside of C + B(0, Vf), C > 0. In the set C + B(0, Vf), the symbol of 
pS,e j g j nc [ e p enc [ en t of e modulo 0{e and is of the form 

P 5 ~ p s + /i 2 r 2 + . . . (2.22) 

where 

~ dist(p,C) 2 , Rep 5 (p) ~ dist(p,C) 2 . (2.23) 
Here we write -6(0, r) ={pG M. 2n ; \p\ < r}, r > 0. 
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As has also been recalled in section 8 in [12], in the set C + B(0, Ve), we have 

while when away from C + B(0, Ve), we shall use that 

p5,e_pS = q (s + p 2 ) . (2.24) 

Here, following [12], we use the notation f £ = 0(e) to express that 

d^f £ (x,o = 0(^-^(0^), 

for arbitrary multi- indices a, f3 G N n . 

We shall study estimates for the real part of the quadratic form associated to the 
operator P s ' £ . The starting point here is Proposition 7.1 of [12]: let < k £ = 0(e) 
be equal to e in C + B(0, y/e) and have its support in C + B(0, \/2e). Let 

K e = Ph (k £ ) 

stand for the Weyl quantization of k £ (x,hO- It is then established in Proposition 
7.1 of 1121 that 



Re ((P 5 ' £ + K £ )u\u) Ch^j \\uf , u G S, 



(2.25) 



when e = Ah, C > is independent of 5, A, and h is small enough depending on 
these 2 parameters. 

Rather than working with the estimate f)2.25p . we shall use that 

Re ((P 5 ' I + K £ )u\u) > - Ch^j \\u\\ 2 , u G S, (2.26) 



which is proved in exactly the same way as in section 7 of [12]. Here we recall that 
e = Ah, e = Ah, A ^> A. In what follows we shall use that the estimate (12.261) holds 
also for u G V(P 5 ' 1 ) = {u G L 2 ; P^u G L 2 } = A~lV(P). 

Let 

n B = nf 



s 



be the spectral projection of P 5 ' e associated with the spectrum of P s,e in the open 
disc D(0,Bh). From Theorem 8.3 in [12] we recall that the spectrum of P s,e in 
-D(0, Bh) is discrete, and the eigenvalues are of the form 

X j>k (h) ~ h (fi jtk + h l ' N i< k n jAl + h 2 ' N ^n j>kt2 + ...), 

where the all numbers in D(0, B) of the form 

for some j G {1, . . . iV}. Here A^, 1 < t < ra, are the eigenvalues of the Hamilton 
map of the quadratic part of p at /3j G C, for which Im A^ > 0. Here we also assume 
that B is chosen such that \(J,j t k\ ^ B, for all j, k. 
Assume that u G L 2 is such that 

MGRan(l-n B ). (2.28) 

We are interested in lower bounds for 

Re(P s '*u\u), (2.29) 

which, in view of (12.261) . amounts to estimating K e u. In doing so, we shall assume, 
for notational simplicity only, that the critical set C defined in (12. lip consists of a 
single point, p\ = (0, 0). From (I2.23p . we know that the leading symbol of P s , p s , is 
such that 

\p 6 (p)\~\p\ 2 , peB(o,Vf). 

Let 

d?d?p s (0,0) 

|a+/3|=2 

be the quadratic approximation of p 5 , so that 

p 5 -p = 0((x, 3 ) = 0({h+ (x, 2 ) 3/2 ) , (s, - (0, 0). (2.30) 

Then p is an elliptic quadratic form on IR 2n , with a positive definite real part. The 
quadratic differential operator 

Po = Op h (p ) : L 2 - L 2 , 
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has discrete spectrum, and from [20] we know that the eigenvalues of Pq are of the 
form hfii t k, with // 1>fe defined as in (I2.27p . 

When estimating K £ u for u G Ran(l — IL3), we also introduce the spectral projection 
H^b associated to Pq and the spectrum of Pq in D(0, Bh). Then, since Ubu = 0, 

K £ u = K £ {U 0iB - U B )u + K e {\ - U 0>B )u. (2.31) 

We shall estimate the first term in the right hand side of (12.311) . using the following 
result. 

Lemma 2.1 We have 

Tl B - n ,B = B (A 3/2 h 1/2 + A- 1 ) : L 2 -> L 2 . (2.32) 



Proof: Let x e Cg°(5(0,2)), < x < 1, be such that = 1 for |(z,OI < 1 - 

Set XVe( x '£) = x(^~ 1//2 ( 2; ) 0)- We shall first show that 



no,,(i-x^) = ^((|) °):^-^ 



(2.33) 



and similarly, that 



n B (l-Xvf)=°B^y)-.L 2 ^L 2 . (2.34) 

When proving (12.331) . we shall use the well-posed Grushin problem for the quadratic 
operator Pq, described in [IT] , [12]. Let 

A .(<^).( 1 + £l±^)!) 1/2 , (,35, 

so that the quadratic elliptic operator Pq is equipped with the natural domain 

V(P ) = {ue L 2 ;A 2 u G L 2 }. 



In section 11 in [TT], using the generalized eigenfunctions of Pq and of the adjoint 
Pq, the authors have constructed the operators 

R . c n _^ L 2 ; R+ ,i} _^ C JV 0j ^ G Nj (2.36) 
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such that when z G -D(0, Bh), the problem 

{P -z)u + R^u_ = v, R + u = v + , (2.37) 

for v G L 2 , v + G C^ has a unique solution it G V(P ), u_ G C^ . Moreover, we 
have the a priori estimate 

h \\A 2 u\\ + \u-\ < 0(1) {\\v\\ + h \v+\) . (2.38) 
Associated to ( 12. 37ft is the Grushin operator 

Vq(z) = (^ P ° r Z R Q y. V{P) x C^ L 2 x C N °, z G D(0, (2.39) 
with an inverse 

fibW = ( ^ ) x _> p( P ) x C^, (2.40) 

depending holomorphically on z. From section 11 of [11], we recall that Sq(z) enjoys 
the following localization properties, when fcGl, 

A 2 ~ k E(z)A k = O : L 2 — > L 2 , (2.41) 

and 

A fc P + (,s) = O fc (l) : C^ -> L 2 , E^(z)A k = O k (l) : L 2 — > C^ , (2.42) 

Let 7 C D(0,B) be a simple positively oriented closed /i-independent contour, 
such that all eigenvalues of Po an d -P 5 ^ in -D(0, Bh) are contained in the interior of 
h'j, so that we have 

dist(z, Spec(Po) U Spec(P^)) > h/0{l), z G hj. 

Here we continue to assume that B > is chosen so that there are no numbers of 
the form fXj^ in (12.27P on the boundary of D(0, B). Writing 

rv = ^ / (z-PoyUz (2.43) 

2tt? J hl 

and using the well-known formula 

(z - Po)' 1 = -E(z) + E + {z)E_ + {z)- l E_{z), (2.44) 
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we obtain that 

n ,B = 7^/ E + {z)E_ + {z)- l E_{z)dz. (2.45) 
2vrz J hy 

Now (12.421) gives that for each fceN, 



E-(z)(l-Xvi)=E_A k A- k (l- Xvf ) = <D U :L 2 ^C^. (2.46) 




Using also that along /ry, we have 

= 0(1) : C N ° -> L 2 , 

and 

J B_ + (^)- 1 = CK/T 1 ) : C^ -> C^ , 

as well as the fact that length of /17 is 03(h), we obtain (12.331) . 

The proof of (12.341) proceeds along the similar lines, relying upon the well-posed 
Grushin problem for P s,£ — z, constructed from the Grushin problem for Pq and 
described in detail in section 11 of [11] and section 8 of [12]. In particular, the 
analogue of the localization property (I2.42p holds true for the inverse of the Grushin 
operator for P S ' £ , and arguing as above, we get ( 12.341) . 

We now come to consider estimates for the difference (IL3 — ^-o,b) X^/fi where we 
claim that 

(U B -U , B )xvf-O B (^ + Vj :L 2 ^L\ e = Ah. (2.47) 

In view of (I2.33H and (12.341) . this will complete the proof of Lemma 2.1. 

When proving (I2.47p . we shall first estimate the difference ^ITb — Ro,b) X^/fi where 

Hb is the spectral projection of the operator P s associated with the spectrum of P 5 
in D(0,Bh). Using I2T43I) together with the similar formula for IL3, we get, by an 
application of the resolvent identity, 

(fi B - n 0iB ) XV f = ^ J (z - P S ) ~' (P 5 - Po) (z - Po)- 1 xvfdz. (2.48) 

Now fl2T4TD . (J232D, (12^41) . together with the fact that along fry we have ||.£r+(z)|| = 
OQi- 1 ), imply that for fcel, 

A fc (z - Po)' 1 A- k = o(^\ : L 2 — > L 2 , (2.49) 
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A 2+k ^ _ p Q yl A -k = Q I 1 ] :L 2^ L 2 ( 2 50 ) 



and even that 

A 2+k (z - Pq)- 1 A~ k = O L 
Here we shall take k = 3 in (12.491) . Writing the integrand in (12.481) as 

(* ~ P 5 ) 1 (P S ~ P°) i* ~ PoV 1 Xvf (2.51) 
= (z - P s ) 1 (P s - P ) A- 3 A 3 (z - Po)- 1 A- 3 A 3 Xvf- 

we see that we have to estimate the operator norm of (P s — Po)A~ 3 . Now it follows 
from (EED, (EUD that 

(P s - P )A- 3 = 0(h 3/2 ) : L 2 -> L 2 . (2.52) 

Also, 

A3 ^=0(^)^ 2 -i 2 - (2-53) 
Combining f[2~3BD . f[2~49]) . (12321) . fl233|) together with the fact that 
( 3 - P 5 )- 1 = O^- 1 ) :L 2 ^L 2 , ze hy, 

which follows from Theorem 8.4 in [T2], and that the length of hy is 03(h), we 
obtain that 

1-J : Z 2 - Z 2 . (2.54) 

It only remains now to estimate the operator norm of ^11 b — Ti-B^j XVf- To that 
end, we write, as in (12.481) . 

(n* - Hb) Xvf = ^-jj z - p5 '") 1 - P ) { z XVfdz (2.55) 

= hi i (* - ( p5,£ ~- ps ) ( z - P T ^ dz 

+ h £ (* ~ P 5 'T (P 5 '"~ P S ) (1 - XV?) (z - Pf 1 Xvfdz 
= I + 11, 
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with the natural definitions of I and II. Using, as in section 8 of [12], that 

(P 5 ^-P S )xvf=o(h^y.L^L\ 

together with the (9(/i _1 )-estimates for the resolvents of P s ' £ and P s along the con- 
tour h'j, we get 

I = O b Qe) = O b (V\ : L 2 - L 2 . (2.56) 
We now come to estimate the term II in (12.551) . We have 



Xvf (P*< £ -P d )=0 lh f ) : L -> L , (2.57) 

and using this estimate, as well as the (9(/i _1 )-resolvent bounds for P 5,e and P s , we 
see that modulo a term whose operator norm on L 2 is 

O b (i)=0 B (A' 1 



we may replace the integrand in II by the following expression, 

(* - P 5 'T X (1 - XVi) (P 5 '*-P 5 ) (1 - Xvf) (z - P 6 )' 1 XV f (2.58) 
Here following (I2.5ip . we write 

(P 5 '"- P 5 ) (1 " Xvf) (z ~ P 5 )' 1 XV? (2-59) 



(p4? _ j») (i _ A-*" 2 (V+ 2 (z - P 5 ) ' A~ fc ) A k Xvf . (2.60) 
Here using (12.241) we see that 

(i*- P) (1 - Xv») A- 1 - 2 = O (^) . (2.61) 



On the other hand, as in (12.501) . 



-i 



A k+2 [z-P d ) A- k = 0[^-):L 2 ^L 2 , (2.62) 



14 



and also, 

A fe X^=0(^J:^->^ 2 . (2.63) 

Combining fl^oTID . fT2ToTj) . (|2^2|) . and (|2^3j) . we see that the expression (jXggp is 
0(1). 

We now come to estimate the remaining factor in (I2.58p . To that end, we let 



h 

be an elliptic symbol in the class defined by the right hand side, and write 
Here we know that 



and since 



(z-P»*)- 1 Lr=o(iy.L>^L> 
LZ' (1 - X J = O (|) : L 2 



it follows that the expression (12.641) is O (=). Using finally that the length of the 
integration contour in (I2.55P is Og(h), we get 

n = B (£) =<D B (l):L 2 ^L 2 . (2.65) 



Combining fl235l) . fl236l) . (j2SS|), we conclude that 



(n fl - n fl ) XV f = O b (V) :L 2 -L a . 



In view of (I2.54p . the bound (12.471) follows, and this completes the proof of Lemma 
2.1. □ 

We now come to estimate the second term in the right hand side of (12.311) . given 
by K £ {1 — Ho,_b)m. Here, the difficulty is that in general, due to a pseudospectral 
phenomenon [3] , the operator norm of Hq,b may exhibit some exponential growth, as 
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h — > 0. To circumvent this issue, our fist task will be to establish a more manageable 
characterization of the vector v = (1 — Hq^u. Specifically, we shall now discuss 
properties of the range of the projection 1 — n 0) _B on L 2 . 

In (I2.27P , following [20J , we have already recalled the form of the eigenvalues of the 
elliptic quadratic operator P . From [20], we know furthermore that the generalized 
eigenf unctions of Pq are of the form 

h -n/4 p (X\ e m*)/h^ (2>66) 



y/hj 

where p(x) is a polynomial on IR n and <$>(x) is a quadratic form with Im$ > 0. The 
degree of the polynomial p(x) in (12.661) tends to oo together with the real part of the 
eigenvalue hfij^ in fl2.27[) . We may also recall from (20] that $ in (12.661) is such that 
the positive Lagrangian subspace A$ = {(x, $'(a;)), x G C n } is the direct sum of the 
generalized eigenspaces of the Hamilton map of po, corresponding to the eigenvalues 
with a positive imaginary part. Correspondingly, the generalized eigenf unctions of 
the formal L 2 adjoint Pq are of the form 

h-»/*q f-^j e l ^/\ (2.67) 

where q is a polynomial and $ is a quadratic form such that Im \I/ is positive definite. 

Let ex, . . . ejsf be a basis for Ran (IIo.b) and let e\, . . . e* N be the corresponding 
dual basis for Ran ((IIo,b)*). If v G L 2 , we have 

N 

Ko,bv = ^(w|e*)e„ (2.68) 

and therefore, v G Ran(l — Ho^b) precisely when v is orthogonal to Ran ((IIo,b)*). 

Proposition 2.2 There exists a selfadjoint h- differential operator Q = Op h (q), 
where q is a positive definite quadratic form on T*~R n , such that 

Ran \Q, ^ C Ran((n , B )*) C R&n(E(Q, CBh)), (2.69) 

for some C > 1 which is independent of Q and B. Here E(Q, A) = 1(-oo,a](Q) is the 
finite rank spectral projection associated to Q and the interval (— oo, A]. 
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Proof: The operator Q will be seen to be essentially the h-Wey\ quantization of the 
classical harmonic oscillator on M n . When constructing Q, recall that the generalized 
eigenf unctions of Pq are of the form (I2.67p . We shall write = (Bx,x), where B 
is a symmetric matrix, B = B\ + iB 2 , where Bj are real, j = 1,2, and B2 > 0. The 
real linear canonical transformation 

k 1 :(x,0^(x,C-B 1 x) (2.70) 

maps the positive Lagrangian subspace A* = {(x, Bx); x G C n } to the positive La- 
grangian subspace {(x, iB 2 x); x G O 1 }. Now since B2 > 0, there exists an invertible 
real n x n matrix C such that the real linear canonical transformation 

^■.(x^O^iC-^C'O (2.71) 

maps {(x,iB 2 x); x G C n } = Ki(A^) to {(x,ix);x G C™}. We take the operator 

3=1 3=1 

associated to = {(x,zx);x G C n } = (tz 2 o ki) (A^), </?(x) = ix 2 /2. To obtain 
the operator Q it only remains to notice that associated to k% and ^2 we have the 
metaplectic operators 

XJ X :L 2 -> L 2 , ^/(x) = e- i(BlX '" )/2 V(x), (2.73) 

and 

U 2 : L 2 —> L 2 , U 2 f{x) = f(Cx) |det C| 1/2 , (2.74) 

both unitary on L 2 , and hence with U := U 2 o C/i, we can take Q := U~ 1 QU = 
OPh ((go (k 2 o Notice that the eigenfunctions of Q are of the form 



&a,h\-E ) &a h=l 



(j^A - e<x,h=i(x) = H a (C 1 x)e i9( - x \ Qe a>h = h (\a\ + ^ e a>h , 

(2 ' 75) 

where H a (x) = YYj=i H aj (xj), a G N n , are the Hermite polynomials. The result 
follows. □ 

Having established a favorable comparison for the linear space Ran ((Ho,_b)*) C 
L 2 , we return to the problem of estimating K £ v, for v — (1 — Uq^b)u. Let ip G 
C°°(R; [0, 1]) with supp(^) C [0, 1], and set 

Mt) = i> (j) , A > 0. (2.76) 
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It follows then from Proposition 2.2 that 

ipBk(Q)v = 0. (2.77) 

To understand the operator occurring in ( 12. 77ft . it is convenient to perform a suitable 
dilation in phase space. Assume therefore that A > in ( 12.761) is such that /i < A < 
1. Let us make the change of variables 

x = A 1/2 x, D x = X~ 1/2 D X . 

Then, since the operator Q is quadratic, 

jQ = ig«(x, hD x ) = Ig« (A 1 / 2 (x, jD^j ^j=q w (x, ^D,)j . (2.78) 

It follows therefore from the functional calculus in the version of [I] that 

i>{\~ l Q) = Op, - [r (x, I £) ) = r w (V 1 / 2 (x, hD x ) ; ^ , (2.79) 

where r G S((-)~ ) for any iV G N, with a complete asymptotic expansion in each 
of these symbol spaces, and with the leading symbol f/>(g(x, £)). 

Remark. It is well known [T3] that when Q = Op™(q) where q is a positive definite 
quadratic form, then the Weyl symbol of f(Q), f G Cg°(R), is of the form f(q; h) = 
f(q) + 0(h). It follows therefore that in ( 12. 79ft we have 

where the leading term of \ is ^{qix, £)). 

It is therefore clear that in (12.251) we can take 

K £ = e?Pbh(Q), e = Ah, (2.80) 

for a suitable choice of B > A fixed, where we can take B as a fixed multiple of A. 
With this choice, we get, using (12.771) . 

K e (i-rv)« = o. (2.8i) 
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Combining fl2T26|) . fl23T|) . Lemma 2.1, and (l2lflT) . we see that for u <E L 2 with 
m G Ran(l — we have 

- Ch^j \\uf < Re (P s > ¥ u, u) + eO B (l 3/2 /i 1/2 + i" 1 ) ||u|| 2 . (2.82) 

Recall that here B > A, B = O(A), is taken fixed, and 5 > is sufficiently small 
but fixed. Choosing first A 3> B large enough and then taking h sufficiently small 
depending on these parameters, we absorb the second term in the right hand side of 
(jM2J) into the left hand side. 

Proposition 2.3 When A < B <^ A, let Us be the spectral projection of P s,e , 
e = Ah, associated with D(0,Bh). Here B is a fixed multiple of A. Assume that 
u G V(P S ' £ ) is such that u G Ran(l — ITb). Then for h sufficiently small, we have 

~ Bh ~ 

Re (P s ' e u\u) > -—r \\u\\ 2 , e = Ah. (2.83) 



Now recall that 



where Ag t e, A s ~ : S — > S, L 2 — > L 2 , have L 2 norm 0^(1). It is therefore clear from 
Proposition 2.3 that if u is such that u G Ran(l — II), where 

n = — [ (z- P)' 1 dz (2.84) 

is the spectral projection of P associated with the spectrum of P in D(0, Bh), then 

\\e- tp/h u\\ < 0(l)e~ t/c \\u\\ , C = C(B)>0. (2.85) 

Therefore, it only remains to consider the restriction of the semigroup e~ tp l h to the 
finite-dimensional subspace Ran(II), generated by the generalized eigenfunctions of 
P corresponding to the eigenvalues of P of modulus < Bh. We shall now proceed to 
do so, in the framework of supersymmetric differential operators. 



3 Supersymmetric operators and return to equi- 
librium in the double well case 

The purpose of this section is to establish Theorem 1.1 in its general form, for a class 
of supersymmetric second order differential operators, including (11. ip . Specifically, 
let 

A : R n -> R n (3.1) 
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be an invertible constant matrix. We decompose 

A = B + C, t B = B, *C=-C, (3.2) 

and assume that 

B > 0. (3.3) 
When 4> e C°°(lR n ; R) is a Morse function such that 

d^(x) = 0(1), <T X ((Bd x <f>, d x <f>)) = 0(1), \a\>2, (3.4) 

we consider the Witten-Hodge Laplacian associated to A and 4>, acting on scalar 
functions, defined as in section 10 of [T2] . 

P = -Af = J2 hD Xj B j>k hD Xk + (M) B i* ^u4>) - hti{B<P") (3.5) 

j,k j,k 

+ $Z c j,khd Xj + hd x . o C iife (0*^)) • 

The principal symbol of P is of the form 

p(x,0 = (Bt,Z) + 2i(C<t>' x ,t) + (B&M, (3.6) 

so that the assumptions ( 12.5}) . ( 12. 6ft . ( 12. 7ft are satisfied. 

Assume that the Morse function has finitely many critical points ii,...ijv el" 
and that 

W(x)\>^, \x\>C. (3.7) 

The assumption (12. 10H holds with 

C = {Pj'J = 1,...N}, p j = (x j ,0), 

and we shall also assume that the dynamical assumptions (I2.14p . fl 2 . 1 5 j) . and (12.161) 
are valid. We then know that the results of section 2 can be applied to P in (13. 5p . 

As in [12], we shall assume now that 

4> has precisely three critical points, of which (3-8) 
two are local minima U±i, and the third one Uq is of index one. 

Then we know from [12] that for C > large enough, P in (13. 5p has precisely 2 
eigenvalues fi = and \i\ in the disc D(0, h/C) for h small enough. Here \L\ is real 
and such that 

\k X = b (a x {h)e~ 2Sxlh + a_i(/i)e" 25 - l/fc ) , Sj = (f)(U ) - <j>{Uj) > 0, j = ±1, (3.9) 
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where aj(h) are real, ctj(h) ~ a^o + ha^i + . . ., a^o > 0. 

The set oo, 0(f/ o ))) has precisely two connected components Dj, j = ±1, 

determined by the condition Uj G -Dj. Let < Xj £ C£°(Dj) be such that Xj = 1 on 
Dj H </> -1 ((— oo, 4>(Uq) — Eq)) for eo > fixed but arbitrarily small. If 

f 3 = /.^(^^(a), J = ±1, 
where Cj(/i) > is a normalization constant such that ||/j|| = 1, and 

f(z-P)- l dz, y=dD(0,^), oo, 
2iri J 7 C 

is the rank 2 spectral projection of P corresponding to the eigenvalues fio = and 
/ii in (13.91) . we have the basis 

e j = n/j, j = ±1, 
for Ran(Il), introduced in [12]. From section 11 of p2] we recall that 

ej = f j + O(h- Nl e-& s '- e0 ')) in L 2 , JVi > 0, 
and that the restriction of P to the space Ran(Il) has the matrix 

( a; 1 )(*->*) = ( Ait' a;a: ) ■ < 310 > 

with respect to the basis (e_i, ei), with the eigenvalues /i = and 

/ii = AliA-i + AiAi. 

A simple computation shows that a corresponding basis of the eigenvectors is given 
by 

Aie_x — A_iei (3-11) 

and 

A^e^+A^ei. (3.12) 

Here we recall from the formulas (11.43), (11.45), and the following discussion in [12] 
that if | Ai| > £7 |A_i| then |A*| > ^ an d AiA^ > 0. We have the same fact 
after permuting the indices —1,1 and the Xj, A*. It follows that 

fii ~ max |Aj| 2 ~ max |A*| 2 . (3.13) 
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Rather than using (13.111) and (I3.12p . we shall make a normalized choice of the 
eigenfunctions, given by 

v = -f= (Aie_i - A_iei) , (3.14) 



lfl X 

and 

vi = —= (AI x e_i + X\ei) . (3.15) 



ffl X 

The corresponding matrix of the coefficients is given by 



V= — [ t ^ )• (3-16) 



1 \ * \ * 



We have det V = 1 and it follows from (13. 13|) that V = 0(1). Hence the inverse 
matrix V~ x has the same properties, so that vq, vi is a well-behaved basis of eigen- 
functions for P. If (e* !, e*) G Ran(n*) is the basis that is dual to (e_i, ei), then the 
corresponding basis of eigenfunctions of P*, dual to (vq,vi) is given by the matrix 
t V~ 1 , so that 

^o* = ( A * e * 1 - A * i e i) > ( 3 - 17 ) 

and 

v*i = -L (A_ i e*_ 1 + A ie *) . (3.18) 
y/Hl 

We summarize the discussion above in the following proposition. 

Proposition 3.1 Let vj and v* be defined as in (13141 . ( 13151) . (13TT7I) . (137131) . Taen 
tae spectral projections 

n, = (»,-, j = o,i 

associated to the eigenvalues fiQ = and /ii m ( 13.91) are uniformly bounded as h — >■ 0. 



Combining (12.851) together with Proposition 3.1, as well as with Theorem 8.4 
of [15] . we get the result in Theorem 1.1 in the general case. 

Theorem 3.2 Let P = -Af where we assume (13. 1113.41) . (13771) . and ~~~) . VFe a/so 
assume that P satisfies the dynamical hypotheses (12.101) . (12.141) . ( 12.151) . so inai tae 
disc -D(0, /i/C) /or C > large enough, contains precisely 2 eigenvalues of P, fi® — 
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and /ii given in (13. 9p . Le£ H, 6e the spectral projection associated with the eigenvalue 
fij, j = 0, 1. Then we have 

Uj = 0(1), /i -»• 0. (3.19) 
VKe /iove furthermore, uniformly as t > ond /i — > 0, 

e -tp/ft = Ho + e -W^ ni + 0{l)e~ t/c , C>0, in C(L 2 ,L 2 ). (3.20) 

Here we have also used that the eigenvalues of P in D(0, Bh)\D(0, h/C) have real 
parts > h/0(l). 

4 Tunnel effect for a well and the sea 

In this section we shall show how to adapt the analysis of section 11 of [12] and that 
of section 3 of the present work to cover the case of a potential with a single well and 
a saddle point, rather than a double well and a saddle point as before. Some parts of 
this section are very close to the corresponding ones of section 8 in [12J , and rather 
than repeating the arguments, we shall often merely refer to the discussion there. 

As in section 3, we shall consider the supersymmetric case. Assume that we 
are given the constant matrices A = B + C and a Morse function 0, that satisfy 
(I3.1l) - (l3.4p . We then have the corresponding Witten-Hodge Laplacian in degree 0, 
given by (13.51) . with a principal symbol (13. 6p so that the assumptions (I2.5p - (I2.7I) 
hold. We refer to the formula (11.3) of [12] for the more general expression for the 
Witten-Hodge Laplacian in degree q > 0, P {q) . 

As before, we shall assume that <fi has finitely many critical points X\, G W 1 
and that I^X 2 -)! — when |x| > C, with C large. The assumptions (12.101) is 
therefore satisfied with C = {pj] j = 1, N} where pj = (xj,0). We also assume 
that the dynamical assumptions (I2.14p . (12.151) . (I2.16P hold. 

As in section 11 of [12J, an application of Theorem 8.3 of [12] to shows that 
the eigenvalues there are of the form 

where jj^ is any eigenvalue of the subprincipal symbol S P ( q ) at (xj,0). From the 
calculations in subsection 10.3 of [12J we recall that the pj t k will be confined to a 
sector {0} U {|argz| <7r/2 — 1/C} around [0, +oo). Recall that it is precisely when 
Xj is of index q (i.e. when the Hessian of at Xj has precisely q negative eigenvalues) 
that one of the pj^ is equal to 0. 
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We shall now introduce more specific conditions for the case that we study here. 
Instead of assuming that we are in the double well case, let us shall suppose that we 
have a single well and a sea, that is 

4> has precisely two critical points, one local ^ 
minimum Ui, and a "saddle point" Uq of index one. 

Notice that this implies that the Maxwellian e~^^ h is no longer an eigenfunction of 
p(°\ since <f>(x) does not go to +00 with 

Put Si = 4>{Uq) — 4>{Ui) so that Si > 0. The set — 00, <j>(Uo)[) has precisely 
two connected components Dj, j = ±1, where Di is determined by the condition 
that Ui G Di, while D_i is unbounded. 

Under these assumptions we shall prove the following result: 

Theorem 4.1 Let P = -A^ be as in ( 13. 5p . where we assume ( 14.21) . Then for 
C > large enough, P has precisely 1 eigenvalue Hi in the disc D(0,h/C) when 
h > is small enough. Here fii is real and of the form 

lii = hai{h)e~ 2Sl/h , (4.3) 

where ai(h) are real, ai(h) ~ 01,0 + + Qi,o > 0, Si = 4>(Uq) — 4>{Ui). 

Remark 4.2. It is clear that Theorem 14.11 implies an analog of Theorem 3.2 in the 
present metastable case. We shall refrain from formulating it explicitly. 

m — — 




level line at 4>{Uq) 

Uo 



Figure 1: A well and the sea 



Proof: We first know that P^ = — Ay* has precisely one eigenvalue [ii = o{h) 
spanning a corresponding 1-dimensional spectral subspace since there is a unique 
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local minimum for 0. Now while e does not belong to L 2 , a truncation of this 
function can be used as a quasimode near U% and it follows therefore as in [12], that 
fii = 0(h°°). Moreover, — A^ has precisely one eigenvalue Jli = o(h) and — Ajf has 
no eigenvalue = o(h) for k > 2 from the discussion in the beginning of the paragraph. 
Since our operators are real we know that the spectra are symmetric around the real 
axis, hence are real. From the intertwining relations 

we then also know that Ji\ = \i\ (see also the discussion at the end of page 69 of 

The construction of the eigenf unctions eo and e\ associated to the critical points 
U\ and Uq for respectively — A^ and — A^ is exactly the same as in [12]. We only 
retain the following from there: 

We begin with -A^ 0) . Let xi e C^pi) be equal to 1 on D x n0 _1 (] ~ °°> 0(^o) - 
eo]) for eo > fixed but arbitrarily small. Consider 

f x = h-^c^tye-^-W^xiix), (4-4) 

where c\{h) ~ ci_o + hc^i + ... > is a normalization constant with cx,o > 0, such 
that || /i|| = 1. Then the normalized eigenfunction associated to /ii is given by 

ei : = f l + C^/r^e"^ 1 -^) in L 2 . (4.5) 

We continue with the study of —A^. Let be the one-dimensional eigenspace 
of corresponding to From an easy extension of [12, Theorem 9.1] (see also 
Remark 9.2 there) to the non-scalar case with the presence of the other non-resonant 
well Ui, we know that E^ 1 ' is generated by an eigenform 

e ( X ; h) = X o{x)e-^ +{x) h-^a {x; h) + 0(e- So/h ), (4.6) 

where Xo £ Co° (neigh (U )) is equal to one near U , S > 0, and 

oo 

a (x; h) ~ y^q 0)fc (x)/i, fc 
o 

is a symbol as in Theorem 9.1 of [12], with ao,o(£^o) 7^ 0. Here the phase 0+ G 
C°° (neigh {Uq)] [0, oo)) satisfies (f)+(x) ^ \x — Uo\ 2 and solves the eikonal equation 
g(x, 0+(x)) = 0, with q = p 2 + p\ - po- 

From [12], let us recall that the phase function 0+ arises as the generating function 
for the stable outgoing manifold through (Uq, 0) for the if^-flow, A$ + , and recall also 
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that </>"([/o) > by Proposition [121 Proposition 8.2]. (Similarly we have a stable 
incoming manifold A^_.) Let k± be the number of eigenvalues of the linearization of 
H q , at that point with ± real part > 0, so that k + + k- = n. Let K + , K_ C 

be the corresponding stable outgoing and incoming submanifolds of dimension k + 
and k_ respectively. Then K + C A^ + , K_ C A^_ and — 0(£/" o ) — <$>± vanishes to 
the second order on ir x (K±). Since 4>"{Uq) has signature (n — 1,1), we conclude that 
dimK + = n — 1, dimfC = 1. It is also clear that A^A^ intersect cleanly along 
K±, so we get 



+ -(0-0(£/ o ))-dist (x,tt x {K + )) 2 , 
- (f>(U ) - 0_ - dist (x, vr x (K_)) 2 . 



(4-7) 



We now make some remarks about the adjoint operator — At^ = (— A^)*" 4 '*. As 
proved in [T2|, Subsection 10.4]), we get the corresponding phases and submanifolds 
in this case, which satisfy 



.* + -(0-0([/ o ))-dist {x,7T x {K* + )f 

0-0([/ o )-0* -dist (sc, 7r x (K* )) 



(4- 



(4.9) 



(4-10) 



and for symmetry reason we recall that 

0- = -0+, 0+ = -0-, 

giving in particular from (14.71) . (14.81) . 

- 0(f/ o ) + 0; - dist (x,^^)) 2 
- (f)(U ) + 0+ - dist (x, ir x (K*_ )) 2 . 

Let /4 be the eigenvalue of Pi " 1 := — that is o(h). As before, this is also 
the eigenvalue o(h) of — A^ and the corresponding eigenspaces Ei 0) and Ei X) are 
respectively spanned by the eigenfunctions 

e*(x; h) = X o(x)e- L ^h-%a* (x; h) + 0{e~ s °' h ) 
and el(x, h) = /r n/ %(/i)e^ w ^ (c/l)) xi(x) + 0(h- Nl e -^ Sl - eo) ) 

Now, using that our eigenvalues and operators are real, we know by duality 
that nl = Hi, and that (Ei°\ E^) and (E» ,E^) are dual pairs for the scalar 
products (u\v) L 2 and {u\v)a respectively. Following [12, Subsection 10.3] we know 
that (aQ (L r o)|a 0i o(^o)) J 4 7^ and that can be normalized so that 

KMa = 1. (4.12) 
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Similarly, denoting by e\ the L 2 normalized eigenfuncion spanning we have 

from fl4.5H4.lin 

(e*| ei ) = l + 0{e-&). (4.13) 

Let (Ax) be the (scalar) matrix of : E^> — > -E'W with respect to the bases (ei) 
and (eo). Let also (A^) be the (scalar valued) matrix of aV* for the same bases. The 
eigenvalue /ii can be viewed as the scalar d^'*d^ : i?( ) — > E^ . We get 

/ii = A*Ai, (4.14) 

and 

Ai = (eSl^ei)x, A* = (fli|d£'*eo)A, j = ±1, (4.15) 

where gi = e*(l + 0(e~ss)) is the vector in E^ that is dual to e\. Here the complex 
conjugate signs are superfluous since we work with real operators, eigenvalues and 
functions. 

We skip the computation of Ai, which is exactly the same as in [12], just recalling 
that the main term is equal to 

- Cl {h)h 1 -* J Xi(x)(A(x)a*(x;/i)|rfx(^))e-^ (0 + (x)+</,(:!;) -* (C/l)) rfx. (4.16) 

and can be evaluated thanks to the stationary phase using (14.101) . We get 

Ai = hk^e-^il + Oie-^)), £i ~£i, + Mi,i + ..., £i,o^O. (4.17) 

similarly, A^ ca be evaluated in a dual point of view as in [I2j and we also get 

A* = hhl(h)e~^ Sl (l + 0(e~3K)), t x {h) ~ f Xfi + ht 1A + t lfi ± 0. (4.18) 

We eventually claim that £i t oE* o > 0. Indeed, this number is real and different 
form zero and if we deform our matrices to reach the selfadjoint case (with A > 0) we 
see that we have a positive sign). Combining this with f)4.14p . the proof of Theorem 
4.1 is complete. □ 



5 Some models of KFP type operators 
5.1 Probabilistic description 

Here we shall give some examples of Kramers-Fokker-Planck type operators. We 
begin with a very short review of stochastic calculus in order to explain their proba- 
bilistic origin, and refer to the books [T7], [T5] for more details. Part of this material 
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can be also found in [TJ, [5] and [6], from where the example of the chain of anhar- 
monic oscillators is taken. 

Let x(t) G M n be a stochastic process satisfying the following stochastic differen- 
tial equation 

dx(t) = b(x(t))dt + adw, (5.1) 

where w is the m-dimensional Wiener process, a is a linear map from M. m to M n , 
and b is a C°°-vector field on M. n , all of whose derivatives are bounded. Under these 
assumptions, there exists a unique global solution x(t) of (15. ip . for a given initial data 
x(0) = x, independent of w, in an adapted stochastic L 2 setting — see the references 
already mentioned. Then we can define a semigroup of operators T', t > 0, by 

R(<p(x(t))\F a )=T t -'<p(x(8)), a.s. (5.2) 

when < s < t. Here T l is the filtration associated to {wis) —w(0); < s < t} 
and x, and tp G C(o)(M n ), where C(o)(M n ) is the Banach space of continuous functions 
vanishing at infinity, with the topology of the uniform convergence. Then T* is a 
strongly continuous positivity preserving contraction semigroup, whose infinitesimal 
generator is given on C^iW 1 ) by 

L = V ■ DV + b(x) ■ V, 

where D = \oa l . The idea now is to extend T* to a larger class of test functions, 
and then to study the evolution of the adjoint (T*)* on the dual space of bounded 
measures. To be precise, let us denote by d\it{x) the probability distribution for x(t), 
defined for allt > 0. We then have 

E(<p(x(t))) = J <p{x)dfi t {x), ^ e C (0) (M n ), 

and we get by ( 15.2ft that /j, t = (T*)*/io, where (T*)* is the adjoint of the operator T* 
acting on the Banach space of bounded measures on M. n . 

We shall now extend the space of test functions. When doing so, we introduce 
the Hilbert space 

H = L 2 (R",e-* (a; W) (5.3) 
where $ G C 00 (R ri ). We shall make the following assumptions concerning $ : 

d^(x) = 0(1), \a\ > 2, (5.4) 
-b(x) ■ V$(x) + -(DV X $, V x $) - -div6 < 0(1) (5.5) 
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and 



d2(h{x)-V^x) + ^(DV x $,V x *)-~tivb) = 0(1), \a\>2. (5.6) 

These conditions will be fulfilled in the case that we shall study in what follows, since 
in the supersymmetric case it is straightforward to verify that they are equivalent to 
(13.41) . Now we can identify the dual TC' of TC with the space of densities 

U* = L 2 {W,e^ x) dx). 

Assume that the measures d\i t are absolutely continuous with density in TC* , and 
write 

dfjtt = f{t, .)dx, 

identifying the measure dfi t with the corresponding density f t . We denote again by 
(T*)* acting on TC* the adjoint of T* acting on TC. We introduce the formal adjoint 
operator L* on TC* of L on TC, with the domain C^°(lR ra ), which is given by 

L* = V • DV - V • b(x). 

We have the following result. 

Lemma 5.1 Assume that $ G C°°(lR n ) satisfies (15. 4p . (15. 6p . Then operator —L on 
7i (resp. —L* on 7i*) is m-accretive, and T l (resp (T 1 )*) is a strongly continuous 
semigroups on 7i (resp. TC*), with infinitesimal generators given by L (resp. L*). 

Proof. It will be more convenient to work in the unweighted space L 2 (IR n ). To 
this end, if G TC, we write (f>(x) = e^^^ipi^x), ip G L 2 . If 

dt<f> = L<p, 

then the equation satisfied by ip is 

M = (e-^ 2 Le^ 2 ) ij, 

so that 

dti> = (d x + d x ®/2) ■ D(d x + d x §/2)4, + b(x) ■ (d x + d x $/2)i/>. 
Let C be a sufficiently large constant. According to (I5.4ti5.6l) . the operator 
- ((d x + d x ®/2) ■ D(d x + d x $/2) + b(x) ■ (d x + d x $/2)) + C 
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has a symbol satisfying the hypotheses (I2.4l) - (l2.7p . Here we may recall that the 
vector field b is bounded on R n together with all of its derivatives. An application of 
Corollary 3.2 in [T2] shows that its maximal closed realization in L 2 coincides with 
the graph closure on S. Coming back to 7i and denoting by L again its maximal 
closed extension, we get that L + C is maximal accretive, and that T* is a strongly 
continuous semigroup, thanks to the Hille-Yosida Theorem. As for the dual semi- 
group on TC*, we also get (for example, using Corollary 10.6 in [IH]), that the same 
occurs for (T*)* and L*. The proof is complete. □ 

From the preceding discussion, we get the equation satisfied by the density / for 
an initial data / G H*, 



where we recall that D = aa t /2. In particular we have dji t = f(t, .)dx = (T*)*/i in 
the space of bounded measures. 

If there exists an invariant probability measure /i^, then its density M. is in 7i*. 
In our present study we shall essentially make the choice C~ l e~® = Ai, but there 
are cases (see e.g. [7], [5]), where it may happen that no invariant measure is known, 
and that another choice of the function $ is necessary. Such a function Ai will be 
called a Maxwellian of the process. Notice that if it exists, it is a O-eigenfunction of 
L* and positive. 

Remark 5.2. Notice that it may also happen that there exists an invariant measure, 
which fails to be finite. We also associated to it a function that we will call again 
Maxwellian (and denote again by M). In that case of course it cannot be normalized. 



Equation (15. 7p is nearly the Kramers-Fokker-Planck type equation that we stud- 
ied in the first part of the paper. In the following sections we shall also do the 
following two things: first we shall exhibit the semiclassical scaling, which corre- 
sponds to the low temperature limit in the models we are going to study later. 
Second, we shall change our unknown by posing / = e~^ x '^ 2 u (forgetting for a 
while the semiclassical scaling), in order to work in the flat space L 2 rather than in 
Ti. Finally, in the three models that we present in the next subsections (Witten, 
Kramers-Fokker-Planck, and the chain of anharmonic oscillators), we shall recognize 
the supersymmetric structure. 








i.e. 




dtf + (-V • DV + V • b)f = 



(5.7) 
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5.2 Witten and Kramers- Fokker- Planck operators 

We begin with the Witten case. It corresponds to an evolution equation with a gra- 
dient field — jWlx) and a diffusion force coming from a heat bath at a temperature 
T. We have 

dx = -jd x Vdt + ^2ifdw. 

Here x G M™ is the spatial variable, the parameter 7 is a friction coefficient, and w is 
an n-dimensional Wiener process of mean and variance 1. With the notation of the 
preceding subsection, we recover an equation of type (15.11) with D = acr t /2 = jTI d 
and b(x) = —jd x V. Equation (15. 7h for the density in this case is then 

dtf - !Tdlf - -yd x {d x Vf) = / dtf - jd x (Td x + d x V)f = 

J\t=o — Jo I J\t=o — Jo- 

Posing T — h/2 and multiplying by h gives the semiclassical equation 

hdtf - ^hd x (hd x + 2d x V)f = 0. (5.9) 

It is then clear that an associated Maxwellian of the process is 

M(x) = e- 2V ^ h . 
Writing / = jv^^u, we obtain from (15.91) that 

hd t u + ^(-hd x + d x V)(hd x + d x V)u = 0. (5.10) 

Here we recognize the Witten operator W = (—hd x + d x V)(hd x + d x V). In the 
notation of (13.ip ~ fl3.5p . it corresponds to a supersymmetric operator with 

A = 1 -I dl <f>(x) = V(x) 

Assumptions of type (13 Ah on V are then fulfilled if 



d a V{x) 



0(1) when \a\ = 2 



0({x) 1 ) when \a\ > 3 



If we also suppose that V is a Morse function with two local minima and a saddle 
point of index one, such that 

|W| > 1/C, for \x\ > C > 0. (5.11) 
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then the dynamical asumptions (I2.14I) - (I2.16I) are satisfied (we skip the proof here, 
which will be given later in the more complex case of the chain of oscillators). In 
particular, in this case we get Theorem 13. 2[ Of course, the corresponding result 
follows also in the case of a single well and the sea, i.e. when V has precisely one 
local minimum and a saddle point (Theorem 14. II and the remark following it). 

We proceed now to discuss the Kramers- Fokker- Planck case, and follow the same 
method. The stochastic equation of type (15. ip comes here from the Newton law 

dx ydt ^p. 



dy = —^ydt — d x V(x)dt + y/2rfTdw 

The parameter 7 is a friction coefficient, and the particle of position x e MJ 1 and 
velocity y G MJ 1 is submitted to an external force field derived from a potential V, 
with w being an n-dimensional Brownian process of mean and variance 1. With 
the notation of the preceding subsection, we therefore have 

D = <to*/2 = ( ° ^ ) and b(x, y) = ( _^ v _ ^ 

The corresponding equation for the density (15. 7p is then 



d t f - iTAyf + d x (yf) + d y (- iy f - d x Vf) = 
f\t=o = fo 

dtf - ld y .(Td y + y)f + ydj - d x Vd y f = 
f\t=o = fo- 

Posing T = h/2 and multiplying by h gives the semiclassical formulation 



(5.13) 



(5.14) 



hd t f - ^hdy.(hd y + 2y)f + yhdj - d x Vhd y f = 0. (5.15) 

A Maxwellian of the process is then 

M{x,y)=C- l e- 2(v ^ +y2 ^/ h 

where C is a normalization constant. If we write / = Ai 1 ^ 2 u, then (15.151) gives 

1 

hd t u + —(—hdy + y) ■ (hd y + y)u + jy ■ hd x u — d x V ■ hd y u = (5.16) 
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This is the semiclassical Kramers-Fokker-Planck operator (11.11) that we studied in 
|12j . and the first part of the present paper. In the notation of section 3, it 
corresponds to a supersymmetric operator with 

A = \(X ^) and <K x >v) = V(x)+y 2 /2. 

The assumptions (13. 4p are fulfilled if 

d a V(x) = 0(1) when \a\ > 2 

As in the Witten case, if we also suppose that V is a Morse function with precisely 
two local minima and a saddle point of index one, and that 

|W| > 1/C, for \x\ > C > 0. (5.17) 

then the dynamical asumptions (I2.14p - (l2.16p are satisfied, and Theorem 13.21 is ap- 
plicable. The corresponding result in the case of a single well and the sea, i.e. when 
V has one local minimum and a saddle point, is also valid (Theorem 14.11 and the 
following remark). 

6 Chains of anharmonic oscillators 

The last example that we give comes from the series of papers [TJ, [5] [B]. It is a 
model describing a chain of two anharmonic oscillators coupled with two heat baths 
at each side. 

The particles are described by their respective position and velocity (xj, yj) G M 2d . 
We suppose that for each oscillator j G {1,2}, the particles are submitted to an 
external force derived from a potential Vj(xj), and that there is a coupling between 
the two oscillators derived from a potential V c (x2 — x{). We denote by V the sum 

V(x) = V^Xi) + V 2 (x 2 ) + V c (x 2 - xi), 

where x = (xi,x 2 ), and we also write y = (yi,y 2 ). By Zj, j G {1,2} we shall denote 
the variables describing the state of the particles in each of the heat baths, and set 
z = (zx, z 2 ). We suppose that the particles in each bath are submitted to a coupling 
with the nearest oscillator, a friction force and a thermal diffusion at temperature 
Tj, (j = 1,2). We denote by Wj, j G {1,2}, two (^-dimensional brownian motions of 
mean and variance 1, and set w = (wi,w 2 ). 
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Figure 2: Oscillators coupled to heat baths 



The fundamental system of equations of type (15. If) is then written as follows, (see 
[7] for more detail concerning the physical constants) 

dx\ = yidt 

dyi = —d xl V(x)dt + z\dt 
dz\ = —"jZ\dt + ~jX\dt — \J2^T\dw\ 
dz2 = —-yzidt + ^x 2 dt — ^2 r yT 2 dw 2 
dy 2 = —d X2 V(x)dt + z 2 dt 
dx 2 = y 2 dt. 



(6.1) 



The parameter 7 is the friction coefficient in the baths. In the (x, y, z) variables, the 
diffusion matrix and the drift appearing in (15.11) are therefore 



D = aa t /2 





7 TJ d 



and b(x) 




(for simplicity we identified T with the 2d times 2d diagonal matrix with coefficients 
Til and T 2 I. The corresponding equation (15. 7p for the density is then 

d t f - lTd 2 J + d x (yf) + d y (-d x Vf + zf) + dMx -z)) = 
ft=o = fo 

where Td\ stands for T x d\ + T 2 d%. We get 





(6.2) 



(6.3) 



d t f - ld z (Td z + (z- x))f + yd J - (d x V - z)d y f 

ft=Q = fo- 

Notice that it is very close to the Witten and Kramers-Fokker-Planck operator. For 
a semiclassical formulation, we pose T\ = ot\h/2 and T 2 = a 2 h/2, and we multiply 
(16.31) by h. This gives 



hd t f + l ai (-hd Zl )(hd Zl + 2{ Zl - x 1 )/a 1 )f 



-hd Z2 ){hd Z2 + 2(z 2 - x 2 )/a 2 )f 



+ (yhd x f-(d x V-z)hd y )f = 0. 



(6.4) 
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At this stage it is difficult to exhibit a Maxwellian. Indeed the existence of an 
invariant measure is a difficult problem solved in some particular case in [7] . Anyway, 
it is clear that the function 

y, z) = V(x) + y 2 /2 + z 2 /2 - zx 

plays a special role, in fact it is the classical energy at temperature 1. We can also 
check that in the case of same temperatures (a\ = a 2 = f a), the function 

is a Maxwellian of the process. We use this function to define the weighted space 
7i = L 2 [e~ 2 ^ I ah dxdydz) as in ( 15. 31) . and in order to work in the flat space L 2 we 
make the change of unknown 

/ = mTu. 

Equation (15.141) becomes 

hd t u + (-hd Zl + -!-0i - xi)^ (^hd Zl + - (z x - xi^j u 
+ J«2 \-hd Z2 + -{z 2 - x 2 ) ) (hd Z2 + ( — - - ] (z 2 - x 2 )) u 



a 

+ (yhd x - (d x V - z)hd y ) u = 0. 
We impose the following condition on the parameter a : 

a > max {01,0:2} /2. 

which corresponds to a semiclassical study at "reference" temperature ah/2 not too 
low. 

Unfortunately we are not able to find any supersymmetric structure in the case 
of different temperatures, since a Maxwellian is not known in this case. From now 
on we therefore stick to the case of identical temperatures T = h/2 so that 

a — ax — a 2 — 1. 

Equation (16. 5p becomes 

hd t u + - {-hd Zl + (zi - xi)) (hd Zl + {z x - x\)) u 



7 
2 

+ (yhd x - (d x V - z)hd y ) u = 



+ \ (~hd Z2 + (z 2 - x 2 )) (hd Z2 + (z 2 - x 2 )) u (6-6) 
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and the Maxwellian was already exhibited M.% = C 1 e 2 */ /l . This equation can be 
written hd y u + Pu = where 

P = J {-hd Zl + (zi - x{)) (hd Zl + (zi - xi)) 

(6.7) 



+ | (~hd Z2 + (z 2 - x 2 )) {hd Z2 + (z 2 - x 2 )) + (yhd x - (d x V - z)hd y ) 



In the notations of Section 3 f l3.1H3.5p . we can write P as a Witten-Hodge laplacian 
P = — A^' with a supersymmetric phase <fi given by 

(f) = § = V(x) + y 2 /2 + z 2 /2 - zx, (6.8) 
and the non-degenerate matrix A — B + C given by 

l ( / d \ t / \ t / 7 d 



A = - \ -I d with 5 = - , C = - \ -I d 

" 1 <yl d / 1 V jl d / ^ V 

















: 


: 










In order to complete the semiclassical study as in [12J, we only need additional 
conditions on the potentials V\, V 2 and V c . It is clear that the conditions 

d a V £ (x) = 0(1) when |a| > 2, with e = 1, 2 and c (6.9) 

imply (13.41) . In view of the definition (16. 8p . it is straightforward that <fi has exactly 
the same number of critical points than V(x) — x 2 /2 with same index. For this it is 
sufficient to notice that there is a natural splitting of the variables for $ given by 

$ = (V(x) - x 2 /2) + y 2 /2 + (z- xf/2. 

We postpone to the end of this section the proof of the following lemma: 

Lemma 6.1 Suppose that V satisfies (16.91) . If in addition V(x) — x 2 /2 is a Morse 
function and there exists C such that 

\dV(x) -x\> l/C when \x\ > C, (6.10) 

then (13. 7p and the dynamical conditions (I2.14H2.161) are fulfilled. 

As a consequence we can apply Theorem 13.21 to operator P: 
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Proposition 6.2 Consider P given by (16. 7p and suppose that V satisfies (16.91) and 
(16.1 Op . Then if the effective potential V(x) — x 2 /2 is of double well type, (two local 
minima and a saddle point of index 1 ), the hypotheses of Theorem 13.21 are fulfilled 
and as a consequence its conclusions apply to operator P. 

Proof. It is straightforward. From the construction of P, hypotheses (13.1113.31) 
are fulfilled. From (16.91) . hypothesis H3.4[) is satisfied. Since the effective potential 
V(x) — x 2 /2 is a morse function of double well type, then (I3.8P is also satisfied since, 
as already noticed, V(x) — x 2 /2 and <f> have the same number of critical point with 
same index. Eventually using Lemma 16.11 and (I6.10p we get that hypothese (13.71) 
and the dynamical conditions (12.14112. 161) are fulfilled. The proof is complete and the 
conclusions of Theorem 13.21 apply to P. □ 

Remark 6.3. Of course the corresponding result follows in the case of a well and 
the sea, ie when V(x) —x 2 /2 has one minimum and a saddle point (Theorem 14. 1 1 and 
the remark after). 

A simple family of such potentials is given for example by the ones for which 
^i(^i) — xf/2of double well type, V 2 (x 2 ) — x\j2 of single well type, and V c sufficiently 
small. Here is an example of such potentials in Id: 



Here x±, x 2 , x' e R. 

Remark 6.4. We did all the computations in the case of 2 oscillators. It is clear that 
the preceding supersymmetric construction works as well in the case of N oscillators 
coupled and with two heat bathes at each side as in [7j. We did not try to verify 
the dynamical conditions in these cases. Eventually recall that the complete study 
at different temperatures seems difficult to treat (see e.g. the recent work by Hairer 
and Mattingly [9] in the case of 3 oscillators). 

Proof of Lemma 16.11 

We only prove here condition (13.71) and (12 . 16[) since the proof of other ones follow 
the same kind of arguments. The proof uses strongly the fact that the symbols 
are with quadratic growth at most. We write p = p 2 + ip\ + po for the symbol of 




1 



cos (a/). 



10 
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the corresponding Hodge Laplacian on 0-forms (minus the constant r yhd/2) where 
denoting (£, 77, () the dual variables of (x, y, z) we have 

T T 
P2 = ^( 2 , Pi = y£ - {d x V - z)r), p = -(z-x) 2 . 

In particular, with the notations of Section 2, we have 

v(x,y,z,d x ,d y ,d z ) = yd x - (d x V - z)d y . (6.11) 

We denote by x\, x^ the critical points of V(x) —x 2 /2, and notice that the critical 
points of $ are (xj, 0, Xj) for j = 1, N. According to definitions (12.111) and (16. lip , 
the critical set C of p is made of the points pj = (xj, 0, Xj, 0, 0, 0) for j = 1, N. We 
also introduce n x (resp. n xyz ) the orthogonal projections on IR^ (resp. Wt%* yz ) from 

Let now e > be a fixed constant. Since V"(a;) — x 2 /2 is a Morse function, it is 
non-degenerate, so that with (16.101) we get the following: there exists C £ such that 

Vx&ir x (C + B(0,e)), \dV(x) - x\ > 1/C £ . (6.12) 

From this result we get easily the following one: there exists C' £ such that 

V(x,y,z)#n Xjytg (C + B(Q,e)), \dV{x) - x\ + \z - x\ + \y\ > 1/C' £ . (6.13) 

This proves (13.71) . 

In order to complete the proof of (12.161) we study in detail the characteristics of 
the flow generated by the vector field is(x, y, z, d x , d y , d z ). We first notice that the flow 
is complete since the gradient of $ is Lipschitz. Let (xo, y$, z ) £ n x ,y,z(C + B(0, e)). 
We denote by (x(t),y(t), z(t)) the integral curves of v for t G R with x(0) = xo, 
2/(0) = yo, z(0) = z and 

x = y 

V = ~(d x V(x) - z) 
z = 0, 

and we have to study 

Po(exp(tu)(x ,yo,z )) = -\x(t) - z(t)\ 2 = -\x(t) -z \ 2 . (6.14) 
We split the study into two cases. Let t G [0, 1]. 
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First case: Suppose \y \ > \d x V(x ) — z \. We write that 

y(t) = yo+ y{s)ds 

Jo 

= Vo- (d x V(x(s)) - z )ds 



(6.15) 

y -t(d x V(x )-z )- I { I V"{x{r))x{r)dr ) ds 



\J0 

t 



y -t(d x V(x )-z )- J (j° V" 



x(r))y(r)dr ] ds 



Since V" is uniformly bounded, and denoting by Cy > 1 a corresponding bound we 
get 



t 2 



sup \y(s) - y \ < t\d x V(x ) - z Q \ + —C v sup \y(s)\. 
se[o,t] 1 se[0,t] 

This implies that on [0, ty] with ty < 1/Cy < 1/Cy 2 we have 



1 t 2 

- sup \y(s) - y \ < t\d x V{x ) - z \ + — C v \y 

1 se[o,t] ^ 

Since by assumption \y§\ > \d x V(x ) — z \ we get 

sup \y(s) - y \ < 2t\y \ + t 2 C v \y Q \ < 3t\y \. 

s£[0,t] 

We can then write that 



x(t) = xo + J x(s)ds = xq + / y(s)ds 
3 

= x + ty + -t 2 \y \(f)(t) 



(6.16) 



with \<f>(t)\ < 1 on [0, iy]. Recalling (I6.14p and using the triangular inequality, we 
get for t e [0, t v ] that if ty < 1/4 

\x(t) - z \ 2 > max ||a; - z \ - 2t\y \, ^\y \ - \x - z \^ (6.17) 

Let now < 6 < ty. We split again the study into two parts: 
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1. If \xq — z \ > 0\y \, we use the first expression in (16.17P and we get that 

\x{t)-z \ 2 > \x -zo\(l-2t/6) > ^\x -z \ on [0,5/4]. 

2. If 9\y \ > \x — z \, we use the second expression in (16.171) and we get that 

\x{t)-z \ 2 >\y Q \{t/2-e)>%\ yo \ on [tv/2M- 



if 6 < t v /8. 

In all cases we get that there exists a constant cy > depending only on Cy such 
that 

\x(t) - z \ 2 > £vmax{|y |, ko - z \, \d x V(x ) - z \} (6.18) 
on an interval of length at least 9. 

Second case: Suppose |y | < \d x V(xo) — z \. As in (" 16. 1 5[) we can write that 

y(t) =y - t(d x V(x ) - zq) - J (J V"(x(r))y(r)dA ds. (6.19) 

Since V" is uniformly bounded, and denoting again by Cy > 1 a corresponding 
bound we get 

t 2 

sup \y(s) -y + s{d x V{x ) -zq)\ < —Cy sup \y(s)\, 



se[o,t] 2 S £[o,t] 



so that 



sup \y(s) -y + s{d x V{x ) - z Q )\ 

s£[0,t] 

t 2 t 2 (6.20) 

< — Cy sup \y(s) -y + s(d x V(x ) - z )\ + —Cy (\y \ +t\d x V(x ) - z \) . 



2 se[o,t] 2 



1/2 

This implies that on [0, ty] with again ty < 1/Cy < l/C v we have 

1 t 2 

- sup \y(s) -y + s(d x V(x ) - z )\ < —Cy (|y | + t\d x V(x ) - z Q \) 

1 se[o,t] z (6.21) 

<Cyt 2 \d x V(x Q )-Z Q l 
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since by assumption \y \ < \d x V(x ) — z \. We therefore get 

sup \y(s) - y + s(d x V(x ) - z )\ < 2C v t 2 \d x V(x ) - z \. 

s£[0,t] 

We can then write that 

x(t) = xq + / x(s)ds = xq + y(s)ds 



o 



(6.22) 



+ ty - ^(d x V(x ) - zq) + ^C v t 3 \d x V(x ) - z \^(t) 



with \ij)(t)\ < 1 on [0, iy]- Recalling (16.1 4p and using the triangular inequality, we 
get for t G [0,t v ] that if t v < 3/(8Cy), 

\x(t) - z \ 2 > max{ \x - z Q \ - t\y \ - t 2 \d x V(x ) - z \, 

t\Vo\ - \ x o ~ zo\ - t 2 \d x V(x ) - Zq\, ^ 23 ^ 

t 2 

— \d x V(x ) - z \ - \x - z Q \ - t\y Q \ } 

Let us now take again < 9 < ty. We consider three cases: 

1. If \xq — z \ > ma,x{9\y \,189 2 \d x V(x ) — z \}, we use the first expression in 
(16.231) and we get that 

\x(t)-z \ 2 >\x -z \(l-t/9-l8t 2 /9 2 )>^\x -z \ on [0,0/6]. 

2. If 9\y \ > max{|a;o — z \, 189 2 \d x V(x ) — z \}, we use the second expression in 
(16.231) and we get that 

\x{t)-z \ 2 >\y \(t-9-t 2 /{lW))> 6 -\y \ on [29,W\. 

3. If 189 2 \d x V (xq) — Zq\ > max(|x — Zo|,0|yo|)> we use the third expression in 
(16.231) and we get that 

\x(t)-z \ 2 > \d x V(x )-z \(]t 2 -18t9-189 2 ) > ^\d x V{x Q )-Zo\ on [t v /2,t v ]- 

4 16 

if 189 2 < t v /16. 

In all cases we get that there exists a positive constant c' v only depending on Cy 
such that 

\x{t) - z \ > c v max {|y Q |, \x - z \, \d x V(x ) - z \} (6.24) 
on an interval of length at least 9. 
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Conclusion From (I6.18p . (16.241) and using (I6.14p . we get that in all cases there exists 
an interval of length at least mm(6,tv/2) on which 

Po (exp(tv)(x ,y ,z )) = -\x(t) - z \ 2 

1 2 

> -(mm(c v ,c v )ma,x{\yo\,\x - z \,\d x V(x ) - z Q \}) • 

(6.25) 

Since (x , yo, z ) ^ K xyz {C + B(0,e)) and using (16.131) we get that the dynamical 
condition ( 12.161) is fulfilled. The proof is complete. □ 
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